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D i Abstract. We consider the gravastar model where the vacuum phase transition 

^0 | between the de Sitter interior and the Schwarzschild or Schwarzschild-de Sitter exterior 

geometries takes place at a single spherical (5-shell. We derive sharp analytic bounds on 
the surface compactness (2m/ r) that follow from the requirement that the dominant 
energy condition (DEC) holds at the shell. In the case of Schwarzschild exterior, the 
highest surface compactness is achieved with the stiff shell in the limit of vanishing 
(dark) energy density in the interior. In the case of Schwarzschild-de Sitter exterior, in 
^q' addition to the gravastar configurations with the shell under surface pressure, gravastar 

configurations with vanishing shell pressure (dust shells), as well as configurations with 
the shell under surface tension, are allowed by the DEC. Respective bounds on the 
surface compactness are derived for all cases. We also consider the speed of sound on 
the shell as derived from the requirement that the shell is stable against the radial 
perturbations. The causality requirement (sound speed not exceeding that of light) 
further restricts the space of allowed gravastar configurations. 
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PACS numbers: 04.40.Dg 
1. Introduction 

The gravastar model has been proposed by Mazur and Mottola (MM) as a possible end 
point of the gravitational collapse of a massive body [U E]. As opposed to the classical 
scenario of the gravitational collapse into a Schwarzschild black hole, the gravastar 
concept avoids the formation of the black hole horizon. Technically, the gravastar 
is a global, static, spherically symmetric solution to the Einstein equations. In the 
interior there is a segment of the de Sitter geometry which is matched to the exterior 
Schwarzschild geometry by means of a spherical shell of matter. In the original MM 
model the quantum phase transition takes place within the boundary layer/shell, but 
the full quantum treatment is immediately replaced by the mean field approximation 
in which (classical) Einstein equations are valid. The de Sitter geometry in the interior 
of the gravastar implies constant positive (dark) energy density p accompanied by the 
isotropic negative pressure p = —p. Such energy-momentum content of the spacetime 
is equivalent to introducing the positive cosmological constant A = 8np into the 
Einstein equations, and can be understood as the gravitational vacuum, or gravitational 
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condensate, which motivated the coining of the term gravastar (gravitational vacuum 
star). In this original model of MM, the phase transition layer consists of two spherical 
5-shells with vanishing surface energy density but with non-vanishing surface pressure, 
between which there is a thin layer of stiff matter (p = p > 0). The transition 
layer is positioned arbitrarily close to where the horizon of the Schwarzschild geometry 
is expected to form. In this way the gravastar model allows for arbitrarily high 
gravitational surface red-shifts and is, for an external observer, indistinguishable from 
being a true black hole. 

A simplified gravastar model was constructed by Visser and Wiltshire [3] in order 
to study the stability properties. The phase transition layer was replaced by a single 
spherical 5-shell. The issue of stability could not be fully settled because the equation 
of state for the shell is not available. However, it was shown that physically reasonable 
equations of state could lead to gravastar configurations that are stable against radial 
perturbations. This simplified gravastar model also motivated a number of authors to 
consider different vacuum geometries in the interior and the exterior. This includes the 
models with anti-de Sitter interior and the Schwarzschild-de Sitter exterior geometries 
and also the Reissner-Nordstrom geometries [Ij, the Born-Infeld phantom gravastar 
[5], and the non-linear electrodynamics gravastar [Sj. A different development of the 
gravastar idea went in the direction of constructing models with continuous profiles for 
the energy density and the anisotropic pressures pD, [8]. In these models the energy 
momentum tensor in the centre of the spherical structure satisfies p = —p and through 
the layer of anisotropic principal pressures (crust) smoothly joins with the exterior 
vacuum spacetime. Gravastars with shells of finite thickness but without 5-shells at 
the boundaries with the vacuum regions have been considered in Ref. [9j and it has 
been pointed out that quasi-normal modes of oscillations might be a possible way of 
distinguishing gravastars from black holes. An account of the observational constraints 
placed on the gravastar model has been given in Ref. [TO] . The gravastar models are 
also closely related to a broader subject of dark energy stars, see eg. Ref. [UJ. In fact, 
the ideas of astrophysical bodies with the de Sitter core representing the gravitationally 
collapsed matter at extremely high density can be traced back to the work of Sakharov 
[12] and Gliner [13] and have been further investigated by Dymnikova in a series of 
papers starting with [TS] . see eg. Refs [T5l [T6] and references therein. 

Given that the physical mechanism behind the proposed vacuum phase transition 
within the gravastar shell is essentially unknown, the energy conditions of GR p2] 
provide the most widely accepted framework of model-of-matter independent constraints 
one can impose on the 'exotic matter' comprising the shell. It is therefore important to 
fully investigate the constraints on the gravastar model that are imposed by the energy 
conditions. In this work we will consider the £-shell gravastar model of Ref [3] and we 
will derive the sharp analytic bounds on the surface compactness that follow from the 
requirement that the dominant energy condition (DEC) is satisfied at the shell. The 
DEC will include the weak (WEC) and the null energy conditions (NEC), but not the 
strong energy condition (SEC). Since the SEC is violated from the outset by the dark 
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energy in the interior of the gravastar, it does not appear natural to require that it holds 
throughout the phase transition layer. 

In Sec. El we review the 5-shell gravastar model and we introduce the configuration 
variables that we use through the rest of the paper. In Sec. [3] we discuss the energy 
conditions of GR imposed on the energy-momentum tensor of the gravastar shell and in 
Sec. H] we carry out the detailed analysis of the gravastar configurations allowed by the 
DEC. In Sec. [5] we discuss the properties of the equation of state as derived from the 
assumed stability of the gravastar with respect to radial perturbations. Conclusions are 
given in Sec. [61 



ds 2 = g a pdx a dxP = -k (1 - p(r)) dt 2 + — dr 2 + r 2 dO 2 (1) 



2. The gravastar model 

The line element for the interior and the exterior geometries can be written using the 
geometrized units G = 1 = c and the coordinates x a = (t, r, <p) in the form 

1 

p(r 

where dVL 2 = d$ 2 + sin 2 $ dip 2 is the metric on the unit 2-sphere, k > is the time 
coordinate scaling constant, and p(r) is the compactness function related to the quasi- 
local mass function m(r) by 

= 2m(r) = 2 f ^ _ 2 ^ ^ 
r r J 

where p(r) is the energy density. In general, to avoid the formation of the event horizon 
within the spherical body the compactness must be less than unity. 

In the interior of the gravastar shell we have the segment of the de Sitter geometry 
with the constant energy density pj nt = Ai nt /87r > and the compactness is given by 

/ \ ^7T 2 -^int 2 „ /o\ 

Pmt{r) = — Pint?" = ~ r > r < a. (3) 

At r = (3/Aint) 1 / 2 the interior compactness reaches unity which corresponds to the 
position of the de Sitter (cosmological) horizon. Therefore the surface of the gravastar 
must be located at the radius r = a < (3/Aint) 1 ^ 2 . 

In the exterior we have the Schwarzschild-de Sitter geometry with the constant 
energy density p cxt = A oxt /87r > and mass parameter M > 0. The exterior 
compactness is given by 

2M 8tt 2 2M A ext 2 

Pext(r) = h — Pcxt r = h — - r , r > a. (4) 

r 3 r 3 

With A ext = 0, M > 0, this is clearly the Schwarzschild geometry, while with A cxt > 0, 

M = 0, this is the de Sitter geometry. In the general case we are considering, A ex t > 0, 

M > 0, the gravastar surface must be located in the region where p ext < 1. For this 

region to exist the condition 9M 2 A ex t < 1 must hold, whereupon the radial coordinates 

of the two horizons of the Schwarzschild-de Sitter geometry which bound the region are 

given by the roots of p cx t( r ) = 1 (a cubic). Thus, if the gravastar surface is located 
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at r = a such that Pi n t( a ) < 1 an d p e xt(oO < 1, no horizon forms within the gravastar 
geometry up to the outer horizon of the Schwarzschild-de Sitter geometry. 

Since the original motivation behind the MM gravastar model was to construct a 
highly compact object that could represent a gravitationally collapsed body of positive 
gravitational mass, one expects the (dark) energy density inside the gravastar to be 
greater than that of the surrounding space, i.e. p int > p ext , or equivalently A int > A cxt . 
We refer to the above as the "gravastar requirement". Note also that the constants 
Ai n t and A ext are introduced for notational convenience; however, A ext can be viewed 
as the cosmological constant A = A ext present in the Einstein equations provided 
that in the exterior of the gravastar shell we set p ext — > p' cxt = and in the interior 
we set pi n t — > p' int = (-^int — A cxt )/87r. Thus, in the picture with the cosmological 
constant, demanding that the energy density in the interior of the gravastar shell is non- 
negative, p' int > 0, is equivalent to our "gravastar requirement". Therefore, although it 
is technically possible to construct the gravastar-like solutions to the Einstein equations 
with p^t < Pext and M > 0, we will not consider them in this analysis. 

Through the rest of the paper we will describe the gravastar configurations in terms 
of the configuration variables defined as follows: 



x and y are the values of the compactness on the interior and the exterior side of the 
shell, while z is the contribution to the compactness at the exterior side of the shell due 
to the energy density in the exterior. 

For the gravastar geometry we require 



z and x are non-negative since the energy densities in the interior and in the exterior 
are non-negative by assumption. The condition z < x reflects the above mentioned 
'gravastar requirement'. The condition y > z ensures that the mass of the gravastar, 
M, is non-negative, while y > x, as will be shown in Sec. |3l ensures that the energy 
density of the gravastar shell is non-negative. 

The gravastar shell is a timelike hypersurface of constant radius r = a along which 
the interior and the exterior geometries are joined. The Israel's thin shell formalism 
[T8] (for textbook coverage see eg. [T9l EDJ) allows one to assign the energy- momentum 
distribution to the hypersurface, as required to make the resulting spacetime a solution 
of Einstein equations. The first junction condition requires that the induced metric 
on the hypersurface be the same for the metrics on both sides of the hypersurface. 
Parametrizing the hypersurface with the coordinates y a = (t, <p), the hypersurface 
metric tensor h a b as induced from the metric (CD) is given by h a b = g^e^e^, where 
e° = dx a /dy a . The hypersurface line element is given by 




y = Pext (a) = 2M/a + A cxt a 2 /3 = 2M/a + z, 
z — A ex ta 2 /3. 



(5) 



0<z<x<y<l. 



(6) 



ds 2 = h ab dy a dy b = -k (1 - p(a)) At 2 + a 2 dft 2 . 



(7) 
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Therefore, for the first junction condition to hold we must have k int (l — /U int (a)) = 
fceoct(l — /w(°))- This is accomplished by setting k ext = 1 for compatibility with 
the usual form of the exterior Schwarzschild metric in case A ex t = 0, and setting 
hnt = (1 — — x ) to appropriately rescale the time coordinate in the interior. The 
second junction condition requires that for the smooth joining of the metrics on the two 
sides of the hypersurface the extrinsic curvature on the hypersurface, K ab = fi a; ^e°e^, 
be the same on both sides of the hypersurface. The nonvanishing components of the 
extrinsic curvature for the hypersurface of radius r = a embedded in the metric of the 
form (EEl) are 



K 



//(a) 



2^1 -M<0 



and 



K* = -v/l-Ma). 

(X 



(8) 



In the gravastar the second junction condition is not satisfied. The joining of the metrics 
is not smooth, but is still possible if one allows for the energy-momentum distribution 
on the hypersurface. The standard expression for the shell energy-momentum tensor on 
a timelike hypersurface is 



st = -^([K]-[KK. 



(9) 



where the brackets denote the discontinuity of the quantity over the hypersurface and 
K = K%. Using the compactness functions for the exterior and interior metric given by 
(J3J) and and the notation (jSJ), the surface energy density can be written 



a 



47r a 



V 1 -M< 



1 



47r a 



X 



V 



(10) 



and the isotropic surface tension can be written 



e = -si 



CP 



57T a 



1 -M°) -r/x'(a)/2 



57T a 



1 



2x 1 



Note that by definition the surface tension has the opposite sign of the surface pressure, 
i.e., 8 > indicates that the surface (shell) is 'under tension', while 9 < indicates that 
it is 'under pressure'. 



3. Energy conditions 

The four standard energy conditions of GR formulated in terms of the components of 
the energy-momentum tensor of the gravastar shell — the surface energy density a ( TTDT) 
and isotropic surface tension 9 (TTTj) — read as follows 

• Weak energy condition (WEC): o > and a — 9 > 0, 

• Null energy condition (NEC): a — 9 > 0, 

• Strong energy condition (SEC): a — 9 > and a — 29 > 0, 

• Dominant energy condition (DEC): a > and a — \9\ > 0. 
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One can see from the structure of (TTOT) and (TTTT) that the status of the energy conditions 
for the gravastar shell can be studied in the space of the three configuration variables, 
x,y,z G [0, 1). As pointed out in Introduction, we require that the energy momentum 
tensor of the gravastar shell satisfies all energy conditions except the SEC. Since the 
DEC is more restrictive than the WEC or the NEC, we proceed to consider the DEC 
through a number of steps. 

We begin with the condition of non-negative surface energy density which, as can 
be easily seen from ( fTUj) . reduces to 

a > x < y. (12) 

This requirement has already been encoded in (J6j). 

Continuing with the condition a — \9\ > 0, we break it into cases according to 9 ^ 0. 
We find 

9^0 =}► x%f (y,z), (13) 



where the function fo is given by 

fo(y,z) = -(l3 + y + Qz- K — } - 

62 \ 1 — y 



y 

- 2 ~ y - 3Z a/36 + y 2 - 6y(6 - z) - 6z(A - 3^)) . (14) 

1 - y ' / 

The configurations with 9 = (vanishing surface tension/pressure), satisfying x = 
fo(y, z), are understood as the shells of dust. The DEC is clearly satisfied for dust shells 
provided (|T2|) holds. For the case 9 < 0, i.e. shells under (positive) surface pressure, the 
condition a — \9\ > reduces to 

9<0, a + 9>0 f (y,z)<x<f+(y,z), (15) 

where /+ is given by 

f+(y, z) = — (61 + 25y + 30z - ^ - ^ 
j + \y, ^ 128 V y l-y 

- 6 - 3Z a/100 + 25y 2 - 9(4 - z)z - 2y(62 - 15^)) . (16) 

For the case 9 > 0, i.e. shell under (positive) tension, cr — |6*| > reduces to 

#>0, a-9>0 =^ f_(y,z)<x<f (y,z), (17) 

where 

f(y,z) = l- , 4(1 xo . (18) 
17 vy ' ; (2-3|/ + 3z) 2 V ; 

Combining the above cases with the condition ( jT2l) allows us to summarize the 
DEC (which includes the WEC and the NEC) with the following conditions on the 
configuration variables x,y,z e [0, 1): 

fo(y, z) < x < f+(y,z) 9 < (pressure) 

-<•<//■ <z = / (y,*) ^ = (dust) (19) 

f-(y, z) < x < f (y, z) 9 > (tension) 
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The configurations saturating the DEC with x = /+ and x = /_ can be called stiff 
and anti-stiff shells. The implications of these conditions on the allowed gravastar 
configurations are studied in the next Section. 

4. Allowed configurations 

The shape of the functions /o, /+ and /_, given by (TUT) . (fl6l) and ( fl8l) . determines the 
regions in the space of configuration variables x,y,z e [0, 1) that are allowed by the DEC 
( |T9l . In addition, we recall the 'gravastar condition' z < x, expressing the requirement 
that the (dark) energy density inside the gravastar is not less than the energy density in 
the exterior spacetime. Before dividing the analysis into cases corresponding to specific 
ranges of the configuration variable z, we observe some general properties of /q, /+ and 
/- 

As already pointed out in Sec. El y < z implies M < (see Eq. (0J), but it can also 
be shown that 

f-(y,z)>f+(y,z), 0<y<z<l, (20) 

i.e., the DEC ( Tl9l) cannot be satisfied with y < z. 

For y = z, which corresponds to M — 0, one can show that 

f-(y, z) = f Q (y, z) = f+(y,z) = z, < z = y < 1, (21) 

indicating that the only y = z configuration satisfying the DEC is x = y = z, which 
further implies Aj nt = A ex t, i-e., simply de Sitter spacetime with no gravastar. 
For y > z it can be shown that the following relations hold, 

fo(y,z) < f+(y,z) <y, 0<z<y<l. (22) 

This indicates that the DEC satisfying configurations with 9 < (shell under surface 
pressure) exist for all z. 

In Fig. [T] the DEC satisfying gravastar configurations are indicated with the shaded 
regions in the x, y plane for several fixed values of z. Details are discussed in what 
follows. 

4.1. z = o 

With z = we are actually considering the gravastar with the Schwarzschild exterior 
geometry. One can show that fo(y,0) < for < y < 1, meaning that in this case 
the dust shells and the shells under surface tension cannot satisfy the condition ([6]). 
Inspection of f + (y,0) reveals that it is greater than zero for < y < 24/25 with a 
maximum at y = 4/5, 

/+(4/5, 0) = ^(19 - Vl05) ~ 0.274. (23) 

It follows that the z = gravastar configurations allowed by the DEC involve the shell 
under surface pressure and are constrained by 

= z < y < 24/25, < x < f + (y, 0). (24) 
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The corresponding region in the x, y plane is shown in Fig. [T] (upper- left plot). 

In the original Mazur-Mottola gravastar model it is assumed that the compactness 
on both sides of the vacuum phase transition layer will approach unity (horizon 
formation) to induce the vacuum phase transition, i.e. x, y — > 1. As we can see from 
the above analysis, the DEC satisfying configurations are clearly pushed away from the 
Mazur-Mottola limit. The highest compactness on the exterior side of the gravastar 
shell that can be reached without violating the DEC is y max = 24/25. It corresponds to 
the surface gravitational redshift Z max = (1 — y m a X ) -1 ^ 2 — 1 = 4, and is obtained with 
flat interior geometry. The same upper bound on the surface compactness was obtained 
in a similar context [21] where the DEC-satisfying shell was considered outside of a 
Schwarzschild black hole of mass m, in the limit where m —>■ 0. The finding which is 
relevant in the context of the gravastar model is that the introduction of the de Sitter 
geometry (dark energy) in the interior, relative to the flat interior geometry (hollow 
shell), does not help to support the DEC-satisfying shell against the collapse at high 
values of the surface compactness, as it might be naively expected due to the repulsive 
character of the de Sitter interior (free falling particle accelerating away from the centre). 
We have also obtained the upper bound on the compactness at the interior side of the 
shell; the highest x that can be reached without violating the DEC is x max ~ 0.274, 
given by (123|) . lowering the maximal exterior compactness to y = 4/5 (surface redshift 
Z = v^5 — 1 ~ 1.236). Let us also note that in the model of Mazur and Mottola the 
vacuum phase transition shell is contributing a negligible fraction to the net gravitational 
mass of the gravastar. Here we see that such a condition cannot be met if the DEC is 
to be satisfied. 



4.2. 0<z<l/3 

With z > we are introducing the Schwarzschild-de Sitter exterior geometry. Let us 
begin with the configurations with the gravastar shell under (positive) surface pressure 
{9 < 0). The functions fo an d /+ tend to — oo as y — > 1, but they have a maximum in 
y with the value greater than z for z < y < 1. The area in the x, y plane bounded by 
f and /+ is the 9 < DEC satisfying region, but we must also observe the 'gravastar 
requirement' x > z. Therefore we look for the intersections of fo and /+ with the x = z 
line and we find 

JoU/o, z) =z, y = yo{z) = _ (25) 

and 

f+(y+, z) =z, y+ = y+(z) = i _ — . (26) 

The DEC satisfying gravastar configurations with the shell under surface pressure 
(9 < 0) are therefore 

, f o < x < f+(y, z) for z < y < yo(z) , x 

0<z<l/3, { -J+Wi ) y-yuv; ^ 

z < x < f + (y, z) for y < y < y + (z) 
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where y and y + are given by (|2"5"|) and (1251) . The dust shells {6 = 0) configurations are 
those with 

< z < 1/3, x = f (y, z), z<y<y Q . (28) 

Upon showing that f_(y,z) < z for z < y < yo(z) it follows that the DEC satisfying 
configurations with the shell under surface tension (6 > 0) are 

< z < 1/3, z < x < f (y, z), z<y<y . (29) 

The DEC satisfying gravastar configurations for the case z — 1/4 are shown in Fig. [1] 
(upper-right plot). 

Considering z > puts the gravastar model into the context where positive 
cosmological constant A is present in the Einstein equations. Here we can no longer 
discuss the surface gravitational redshift, but we still find it important to discuss the 
maximal surface compactness that can be reached without violating the DEC. As in the 
z = case, the highest compactness on the exterior side of the shell is achieved with the 
'hollow interior', x = z, which here means that the (dark) energy density in the interior 
is equal to that in the exterior; the maximal compactness is y+(z) given in (1261) . It can 
be shown that y+(z) starts with the value 24/25 at z = and increases toward unity 
as z — ► 1/3. However, if the (dark) energy density in the interior is higher than that in 
the exterior region, x > y, the maximal compactness that can be achieved is lower than 
that given by y + . Another interesting feature of gravastar configurations with z > 
is that in addition to the configurations with shell under surface pressure here we can 
have shells of dust and shells under surface tension. 

4.3. z = l/3 

The z — 1/3 case can be seen as the critical case between the regimes z ^ 1/3. The 
functions fo(y, 1/3) and f+(y, 1/3), bounding the 9 < DEC satisfying region, are 
monotonically increasing in y for 1/3 < y < 1, reaching finite values /o(l, 1/3) = 1/2 
and / + (1, 1/3) = 3/4, while f-(y, 1/3) is monotonically decreasing diverging to — oo 
as y — *■ 1. Observing also the 'gravastar requirement' x > z, the DEC satisfying 
configurations are simply 

(Mv,z)<x<f + (y,z) (0<O) 
1/3 = * < J/ <1, lx = f {y,z) (6 = 0) (30) 

(z<x<f (y,z) (0>O) 

and are shown in Fig. [1] (lower-left plot). The distinctive feature of this special case 
z = 1/3 is that there is a wide range of configurations that allow the surface compactness 
on the exterior side to approach unity arbitrarily close. The compactness on the interior 
side of the shell is still bounded from above by f + , or fixed by / for a dust shell. For 
example, the dust shell with y — ► 1 will have x — > 1/2, which implies interior (dark) 
energy density which higher by factor x/z — > 3/2 than the (dark) energy density in 
the exterior. However, the value of z = 1/3 considered here implies that the radius 
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of the gravastar is larger than one half of the radius of the exterior horizon of the 
Schwarzschild-de Sitter geometry, which makes these configurations highly unlikely to 
be relevant in the astrophysical context. 

44. 1/3<z<1 

For completeness here we consider also the 1/3 < z < 1 configurations. The functions 
/_, fo and /+ tend to 1 as y — > 1. Since /o > z for 1/3 < y < 1, the configurations with 
dust shells and shells under surface pressure (6 < 0) are 

l/3<z<y<l, f {y,z)<x<f+{y,z). (31) 

The shell under tension (8 > 0) configurations must be treated more carefully since it 
can be shown that the function f-(y, z) dips into the x < z region for 1/3 < z < 2/3, 
which is forbidden by the 'gravastar requirement' z < x. We find 

f-(y-, z) =z, y- = y-{z) = _ - , (32) 

so for 1/3 < z < 2/3, the 9 > DEC satisfying gravastar configurations are 

l/3<,<2/3, { * ^ £ JV > (33) 

where ?/_ is given by fl32|) . For 2/3 < z < 1 the 6 > DEC satisfying gravastar 
configurations are 

2/3<*<3/<l, f-(y,z)<x<f (y,z). (34) 

The special case z = 1/2 is shown in Fig. [1] (lower-right plot). The interesting feature of 
z > 1/3 configurations is that the compactness on both sides of the shell can approach 
unity without violating the DEC. However, as already commented in the case of z = 1/3 
configurations, the size of such object is too large to represent a body of astrophysical 
interest. 



5. A note on the equation of state 

In a perfect fluid described by an equation of state of the form p = p(p), where p is 
the isotropic pressure and p is the energy density, the speed of sound propagation, c s , is 
given by the relation q? = dp/dp. In order for the causality to be preserved, it is natural 
to require that the sound speed does not exceed the speed of light, i.e. in our units 
c s < 1. For a view on the status of the relation among causality and the speed of sound 
see Ref. [22]. Here we will investigate the speed of sound propagating (tangentially) 
along the gravastar shell, assuming that the equation of state for matter comprising the 
shell can be given in the simple (barotropic) form 9 = 6(a). Note that while the matter 
comprising the shell can be considered as the perfect fluid on the hypersurface, in the 
complete spacetime it represents a highly anisotropic structure (the transverse pressure 
is described by the ^-distribution) . 
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Lz = 1/4J 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



x x 

Figure 1. Gravastar configurations allowed by the DEC shown as shaded regions in 
the x, y plane for fixed values of z = 0, 1/4, 1/3, 1/2. (x, y, z are defined in ((5])) 

In Ref. [3] a procedure has been developed that allows extracting of the equation 
of state 6 = 0(a) of the matter comprising the gravastar shell on the basis of the 
requirement that the shell be stable against radial perturbations. From the dynamical 
version of the expressions for the surface energy density and the surface tension of the 
shell, which involve the terms a 2 and a, overdot indicating the derivative with respect 
to the proper time of the observer co- moving with the shell, the 'classical' equation of 
motion of the form |d 2 + V = E follows. The 'potential' V is a function of the gravastar 
radius a and other configuration variables, while the 'energy' E = (one is not allowed 
to adjust the energy as in the classical context). Upon identifying d 2 = —2V(a) and 
a = —V'(a) (See Ref. [3] for the full derivation) the dynamical expressions for o and 9 
become functions of the dynamical gravastar radius. The requirement for the stability 
of a shell at the radius a can then be expressed as the requirement imposed on the 
potential, 

Via) = 0, V\a) = 0, V"{a) > 0. (35) 
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Therefore, if one chooses a particular shape of the potential V(a) leading to stable 
gravastar shells, the equation of state 9 = 9(a) can be parametrically extracted by 
varying a over the range of interest. The important special case is V(a) = which, as 
argued in Ref. [3j, reflects the notion of strict stability of the shell. If the shell would 
be displaced to a nearby radius, it would find itself there in a new state of equilibrium. 
In the case V(a) = the dynamical equations for surface energy density and surface 
pressure of the shell coincide with (reduce back to) the static expressions for a and 9 
given by ( flOi) and ( {Til . Therefore we can use this approach to compute the quantity 



= _ d o = _d*/d« (36) 

s do- da/da 



directly from ( flOl) and ( fill) . 

We now proceed to compute c 2 for the gravastar models we have considered so far. 
Using (|36|) . (|T0|) and (TTTT) one can obtain an expression for c 2 in terms of a, A int , A cxt 
and M which is rather involved, but with the use of our configuration variables x, y, z 
defined in ([5]) it becomes much more compact and can be written 

2 _ Av 3 - u 3 [l -6z + 3(v 2 + z)(v 2 + 3z)} 

° s ~ 4u 2 v 2 [-2v + u{-l + 3v 2 + 3z)] ( } 

where u 2 = 1 — x and v 2 = 1 — y. The quantity c 2 diverges for u — > and v — > 
and also at u — 2v/(3v 2 + 3z — 1). To obtain the boundaries of the regions in the 
x, y plane (for the given values of z) in which < c 2 < 1 one must solve c 2 = 0, 1 
for u which leads to a cubic equation. The resulting expressions are very complicated 
so we only plot the results for z = (Schwarzschild exterior geometry) and z — 1/3 
(case considered in Sec. I4.3[) in Fig. [2j The region where causality is preserved with the 
equation of state having usual properties (surface pressure increases as energy density 
increases) are indicated by 'A'. We can see that the causality requirement combined 
with the energy conditions considered before imposes still stronger upper bounds on the 
gravastar exterior surface compactness y. As an example, for the x = z = case (hollow 
shell, Schwarzschild exterior) the highest allowed surface compactness is y max — 0.848 or 
analytically 2(106 + a - 89/a)/255 where a = (225v / 106 - 2195) 1/3 . In both examples 
shown in Fig. [2] we also see that the region 'C where c 2 is negative includes some 
regions that are allowed by the energy conditions. Negative c 2 means that the surface 
pressure decreases as the surface energy density increases, which indicates instability 
with respect to perturbations. Therefore, the requirement that the dependence of 9 
on a is such that stability with respect to radial perturbations is obtained can lead 
to configurations where the gravastar shell is unstable against tangential perturbations 
such as sound propagation. 

6. Conclusions 

Compact objects like gravastars envisaged as possible alternatives to black holes 
are certainly not less controversial than black holes themselves |23j . The essential 
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Figure 2. Behaviour of the quantity Cg = — d8/da (speed of sound squared) as derived 
from the requirement for the stability against the radial perturbations for the gravastar 
shells with z = (left plot) and 2 = 1/3 (right plot), shown in the x, y plane. Dashed 
lines indicate configurations where — > ±oo, while the solid lines bound the region 
'A' in which < cj < 1. Region 'B' is the 1 < Cg superluminal regime, while 'C is 
the Cg < unstable regime. Shaded regions are the gravastar configurations allowed 
by the DEC (See Fig. [J). 

phenomenological feature of a gravastar that makes it difficult to distinguish such a 
body from a true black hole is its high surface redshift. Also, the analysis of radiation 
from a supermassive object Sgr A* [21] implied the existence of a horizon in this object. 
However, the absence of detectable heating could also be consistent with a gravastar 
possesing sufficiently large heat capacity [10J. Another feature which may help in this 
respect is the analysis of quasi-normal modes of oscillation of a gravastar. They are 
found to be different from those of a black hole [9], and could, in principle, provide a 
way to observationally discern a gravastar from a black hole. In the same spirit it is 
important to analyze the constraints on possible realisations of gravastar models. Since 
the physical theory behind the quantum vacuum phase transitions which are assumed to 
take place within the gravastar shell has not been fully formulated, the energy conditions 
of GR provide the most general model-independent constraint on the structure of the 
energy-momentum tensor of the gravastar. We have carried out a detailed analysis of 
the constraints that the energy conditions impose on the possible configurations of the 
single 5-shell version of the gravastar model [3]. 

The most relevant constraint is the upper bound on the surface compactness 
fj, = 2m /r of the gravastar which was in the original MM model assumed to approach 
unity arbitrarily close. Exactly this feature gives raise to arbitrarily large surface 
gravitational redshift which makes the gravastar observationally difficult to distinguish 
from the true Schwarzschild black hole. (Note that the original MM model violates the 
DEC at its two 5-shells that are assumed to possess nonvanishing surface pressure but 
vanishing surface energy density. For a general discussion of known violations of energy 
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conditions see eg. Refs [T§| 125])- 

In this paper we have derived the bounds on the surface compactness that follow 
from the requirement that the matter comprising the gravastar shell satisfies the 
DEC. The procedure is fully analytic and has been made possible by the use of the 
configuration variables x, y, z based the values of the compactness at the inner and the 
outer side of the shell defined in (jSJ). The space of gravastar configurations satisfying 
the DEC is confined within the range f_(y,z) < x < f + (y,z), where the functions /_ 
and / + identifying the configurations with anti-stiff and stiff shells are given by ( 1T81) and 
(|T6|) . In case of the Schwarzschild exterior geometry the maximal surface compactness 
/•fmax = Umax = 24/25 is achieved with a stiff shell and flat geometry in the gravastar 
interior (compactness in the interior side of the shell x — 0), i.e. in the limit of vanishing 
(dark) energy density in the interior. Introducing the (dark) energy density in the 
interior (x > 0) does not lead to higher surface compactness that can be reached. This 
finding may be seen as counter-intuitive since one might have expected that the repulsive 
character of the de Sitter geometry will help to support the shell against the collapse, 
however, as our analysis has shown, this is not the case if the shell is to satisfy the DEC. 

For completeness, we have extended our analysis to the case of the Schwarzschild- 
de Sitter exterior geometry which places the gravastar model in the context where the 
positive cosmological constant A is present in the Einstein equations. Relevant bounds 
on the surface compactness have been obtained for all values of the configuration variable 
z = Aa 2 /3, a being the gravastar radius. For z < 1/3 the highest surface compactness 
is bounded below unity and is achieved with the stiff shell. At z = 1/3 there is a wide 
range of configurations where the surface compactness can approach unity arbitrarily 
close and it includes shells under surface pressure, dust shells and shells under surface 
tension. For z > 1/3 the highest surface compactness is again bounded below unity and 
is achieved with anti-stiff shells. 

We have also applied the procedure of Ref. [3] to compute the quantity c 2 = — AO /da 
which in the context of perfect fluids corresponds to the speed of sound. In order not to 
violate the principle of causality one normally requires that the sound speed does not 
exceed the speed of light. This requirement has given still stronger upper bounds on 
the gravastar surface compactness as compared to the bounds derived from the energy 
conditions. However, we must emphasize that the procedure is based on the assumption 
that the exotic matter comprising the gravastar shell obeys a barotropic equation of 
state of the form 6 = 6(o~), on the assumption that the shell is stable against the radial 
perturbations, and most importantly on the assumption that the fluctuations in surface 
tension and surface energy density propagate along the shell at the speed given by the 
usual formula stated above. Since the full theory describing the gravastar shell is not 
available none of these assumptions can be taken for granted. Similar situation has been 
encountered in the context of the stability analysis of wormhole solutions in Ref. [26] . 

All the above results were derived using the 5-shell gravastar model, but they may 
be relevant in the broader context of the dark energy stars, or gravastars with continuous 
pressure profiles, where one attempts to construct astrophysically plausible models of 
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self-gravitating objects with dark-energy cores. As shown in Ref. [TJ, and applied in 
in Ref. [Sj, such objects require anisotropic principal pressures. The analysis of the 
maximal surface compactness of anisotropic spheres is considerably more complicated 
than in the isotropic case; probably the most general results obtained in this area is 
Ref. |27j . although still not general enough to include the case of the dark energy interior, 
or Ref. [28J where the influence of the cosmological constant on the compactness and 
other physical parameters of anisotropic configurations was considered. However, thin 
shells (soap bubbles) appear to be the most efficient configurations in attempts to reach 
high surface compactness of anisotropic spherical bodies. The maximal compactness of 
the 'shell around a black hole' of Ref. [21J, or that of our flat interior gravastar (hollow 
shell), clearly saturates the general bound derived in Ref. [27], and similar suggestions 
have also been given in Ref. [29] . 
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